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Measurements were made at Hanover, N.H., of the effect of the solar eclipse of 20 July 1963 
on the phase and amplitude of two VLF signals propagated over relatively short paths. The skywave 
signal from NAA, Cutler, Maine, (14.7 kc/s) was reflected once from the D region at a point of total 
optical obscuration. The skywave signal from NSS, Annapolis, Md., (22.3 kc/s) was reflected from a 
point where the maximum optical obscuration was 88 percent. Changes in the phase and amplitude 
of the signals during the eclipse ranged from about 40 percent to 100 percent of the diurnal variation. 
The times of the maximum phase and amplitude effects were used to calculate effective recombina- 
tion coefficients for the ionosphere. The maximum amplitude of VLF signals occurred before the 
time of maximum optical obscuration. This can be attributed to a specific source of ionizing radiation 
above the sun's west limb, which would be uncovered before the optical maximum. The calculated 
effective recombination coefficient is about 2 X 10~ 5 . This value is closer to the value of approximately 
3 X 10~ 6 , calculated from theory by Crain [1961] for the nighttime ionosphere, than to Crain's daytime 
value. Most values of recombination coefficients from previous radio experiments are closer to 
Crain's daytime values. 



1. Introduction 

The total solar eclipse of 20 July 1963 provided a 
rare opportunity for a measurement of the effect of 
the eclipse on a VLF signal. The phase and amplitude 
of a VLF signal propagating over a short path are very 
sensitive to changes in the D region of the ionosphere, 
such as those caused by an eclipse; the sensitivity is 
considerably less for long paths. 

The changes in the phase and amplitude of a VLF 
signal during an eclipse may be used for several 
purposes. First, the changes in phase and amplitude 
may be used to determine the changes in the equivalent 
phase height of reflection and the changes in the 
apparent magnitude of the reflection coefficient. 
The usefulness of these quantities is limited however, 
since it is difficult to relate them to anything other than 
the normal diurnal variations. However, the 
changes in the amplitude and phase may be used to 
determine a reaction time and perhaps an effective 
recombination coefficient from parts of the D region 
which affect the signals. These are important quanti- 
ties, since so little is known about the ionization 
production and removal processes in the D region 
[Reid, 1964]. 

The path of total eclipse in the D region (60-100 
km height) passed about 200 km southwest of NAA, the 
high power U.S. Navy VLF transmitter at Cutler, 
Maine, which was at that time transmitting frequency- 
stabilized signals at 14.7 kc/s. A receiver was set 
up at the Dartmouth College field site near Hanover, 
N.H. The strongest skywave signal which arrived at 
this receiver from NAA came from a single ionospheric 



reflection in the area of the D region in which the sun 
was totally eclipsed. 

In addition, the 22.3 kc/s signal from NSS at 
Annapolis, Md., was recorded. In this case the sun 
was only 88 percent eclipsed at the reflection point. 
Some characteristics of these paths are shown in 
table 1. An outline map of the paths is shown in 
figure 1. 

Table 1. VLF signals recorded at Hanover, A/.//., during the eclipse 
of 20 July 1963 



Transmitter 
call, freq. 


Location 


Path 
length 


Maximum 
eclipse % mag. 


Type of 
measurement 


kc/s 
NAA 14.7 

NSS 22.3 


Cutler, Maine 
Annapolis, Md 


km 
406 

633 


100 
88 


Abnormal 
loop a 

Normal loop a 



a See section on experimental procedure. 

2. Equipment 

The equipment used in this experiment was de- 
signed to record the phase of VLF Standard Fre- 
quency Transmissions relative to the phase of the 
output of a rubidium vapor frequency standard, and 
the amplitude relative to a calibration voltage. The 
antennas were multiturn shielded loops, 1 m square, 
tuned to the received frequency. The system output 
was a pair of d-c voltages for each channel, proportional 
to the amplitude and phase of the coherent signal. 
Each pair was recorded on a dual-channel chart 
recorder. 
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NAA 14.7 kc/s 
CUTLER, ME. 




Figure 1. Outline map showing the VLF paths, the time and width 
of optical totality in the ionosphere, and the times of optical 
maximum at the reflection points. 



3. Experimental Procedure 

It is generally a very good assumption that a VLF 
transmitter produces only parallel polarized 1 signals. 
However, some of the parallel polarized energy is 
converted to perpendicular polarized energy upon 
reflection from the ionosphere, because of the action 
of the earth's magnetic field. This reflection may be 
expressed by the use of the two reflection coefficients 
1(7? 1 1 and _j |/? '■_(_. The two subscripts indicate the polari- 
zation of the incident and reflected ray, respectively 
[Budden, 1951]. The amplitude of the signal reflected 
in the parallel polarization is generally several times 
the amplitude of the converted signal. A normal 
loop has maximum sensitivity to the parallel polarized 
groundwave and skywave signals. An abnormal loop 
has maximum sensitivity to the abnormal skywave 
[Best et al., 1936]. 

The ground distance from NAA to the receiver was 
about 406 km. At this distance, the groundwave is 
still of greater amplitude than the skywave signal. 
Since the groundwave is independent of the eclipse 
effects, it is advantageous to reduce the amount of 
the groundwave present in the recorded signal. This 



is often done by taking advantage of the properties 
of the abnormal loop. The groundwave is initially 
reduced by rotating the loop to an azimuth perpendicu- 
lar to the azimuth of the transmitter. A more precise 
adjustment may then be made by utilizing the fact 
that the ionospheric signal is minimum at local solar 
noon [Best et al., 1936] . The loop may therefore be 
rotated to the point where the signal is a minimum at 
this time to obtain a position close to the optimum for 
suppression of the normally polarized signals. 

The usual procedure is then to observe the diurnal 
variation in phase and amplitude and to make incre- 
mental changes in the position to determine the 
optimum orientation. However, it was not possible 
to observe the diurnal variations for more than a day 
before, and one or two days after the eclipse. The 
position obtained by rotating the loop to a minimum at 
local noon therefore had to be used. Considerable 
processing of the data is necessary to remove the 
remaining groundwave component. 

4. Measurements 

Photographs of the original records obtained during 
the eclipse are shown in figure 2 following a theoretical 
curve showing the fraction of the surface area of the 
sun that is visible at the one-hop reflection point in 
the D region for the NAA-Hanover path. Note that 
the direction of amplitude increase is downward. 

The unobscured fraction of the solar disk was also 
calculated for the extreme points at the ends of the 
major axes of the elliptically shaped first Fresnel zone 
at heights of 60-100 km for both the NAA and NSS 
signals. The results, shown in table 2, indicate the 
spread in magnitude of maximum obscuration and 
times of its occurrence, which would result in an ob- 
scuration curve broader than those shown in figure 2 




NSS (ANNATOLIS.MW-HANOVER^^SO kc/s 




360° 
270° 
180° 
^90° 



ir.Pli 



1 The polarization of the signal is called parallel or "normal" when the electric vector FIGURE 2. Unobscured fraction of Solar disk for NAA reflection 

is in the plane of incidence, and perpendicular or ''abnormal" when the electric vector is • rf y LF records taken at Hanover, N.H., during the eclipse 

perpendicular to it. The terms normal and abnormal loop refer to a vertical loop r on J 1 m^Q 

antenna oriented with its plane in, or perpendicular to, the plane of incidence of the signal. Of 20 July 190d. 
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for the reflection point. The table shows that, for 
the NAA reflection area (first Fresnel zone), the times 
of maximum optical obscuration are spread out over a 
period of about 2 min, and the unobscured fraction is 
0.0142 or less during this period. 

The maximum obscuration for the NSS reflection 
area is spread over about 4 min (—3 to -hi min) 
around the maximum at 2149 UT at the reflection 
point. The values of the unobscured fraction deviate 
up to ±0.04 from the value of 0.12 at the reflection 
point. 

Some characteristics of the equipment should be 
considered when examining the phase and amplitude 
records shown in figure 2. First, because both stations 
were transmitting Morse code traffic, the amplitude 
was averaged with a time constant of the order of a 
second or more. The recorded amplitude level of such 
a signal therefore depends on the average keying rate 
of the signal. The VLF stations transmitted a short 
pulse every few seconds for a period of about 3 min 
in each hour. This is evident on the NAA record in 
figure 2 at about 49 min after each hour. The average 
level drops almost to zero during these periods. 

Immediately following these pulsing periods, there 
is a period of continuous (locked-key) transmission; 
following this is a period of no transmission. The 
recorded amplitude reaches a maximum and a mini- 
mum, respectively, at these times. 

The phase scale is shown at the right of each 
record. When the changing phase reaches zero or 
360° at the edge of a chart, it then moves over to the 
other side of the chart and continues in the same 
direction. The phase trace generally wanders during 
the period of no transmission due to noise or a slight 
unbalance in the system. 

It will be noted that the amplitude recording of NAA 
went off scale during the eclipse. A series of readings 
of the amplitude were taken independently with other 
instruments and these values were used where the 
chart reading was off scale. 

Although, unfortunately, NSS was off the air at 
about the time of maximum effect of the eclipse, a pair 
of smooth curves were drawn through the available 
data points for the amplitude and phase, which are 
believed to represent the actual variation fairly well. 
The fact that the results for phase and amplitude on 
both signals are quite consistent indicates that the 
deduced curves are fairly accurate. 



4.1. Removal of Groundwave Components 

Polar plots (figs. 3 through 8) were made from the 
chart recordings after the drift due to the small dif- 
ference in frequency between the transmitting and the 
receiving reference oscillator had been removed. 
A phase advance is counterclockwise in all cases. 
Relative amplitude from the record origin is shown in 
units of chart deflection that would be attained for a 
continuous signal. 

In order to obtain the variation of the one-hop 
skywave signal, it is necessary to correct for ground- 
wave signal present. This may be done by finding a 
new origin for the polar plot, also shown in the figures, 
which is close to the origin of the skywave vector. 
Two methods were used to determine this new origin. 
The first involved overlaying all the polar plots for a 
given signal and determining an area within which the 
skywave origin probably lies. The curvature, the 
maximum and minimum points, and the overall sym- 
metry of the plot must be considered when using this 
method. 

At the distances involved (400-650 km) the polar 
plots are complicated by a diurnal phase shift of less 
than 360°, rapid fluctuations of phase and amplitude 




FIGURE 3. Polar plot of phase and amplitude of IS A A (14.7 kc/s) 
received at Hanover, N.H., 21 July 1963. 

UT is shown on the curve. The scales are drawn to the origin of the original record. 



TABLE 2. Conditions at optical eclipse maximum at heights of 60-100 km for one-hop reflection 
point and extreme points of first Fresnel zone for paths to Hanover 



Signal 


Reflecti 


on point 




Extreme points of first Fresnel zone 




























UTof 
max. 


% Unob- 
scured 


NE point 
•UT of max. 


% Unob- 
scured 


SW point 
UT of max. 


% Unob- 
scured 


NW point 
UT of max. 


% Unob- 
scured 


SE point 
UT of max. 


% Unob- 
scured 


NAA 
14.7 kc/s 


2145 


0.0 


2145 


0.00 
-.36 


2145 


1.42 « 
-0.66 


2144 


0.00 
-.40 


2146 


0.40 
-.04 


NSS 
22.3 kc/s 


2149 


11.8 
-10.7 


2146 


6.8 
-6.0 


2150 
(60-95 km) 


16.3 
-15.4 


2148 


11.5 
-10.5 


2149 


11.5 
-10.5 
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Figure 4. Polar plot of phase and amplitude of NAA (14.7 kc/s) 
received at Hanover, N.H., 19 July 1963. 

UT is shown on the curve. The scales are drawn to the origin of the original record. 




90° 



Figure 5. Polar plot of phase and amplitude of NAA (14.7 kc/s) 
received at Hanover, N.H., 20 July 1963. 

UT is shown on the curve. The scales are drawn to the origin of the original record. 




2 3 4 5 6 7 8 9 
UNITS OF CHART DEFLECTION 



Figure 7. Polar plot of phase and amplitude of NSS (22.3 kc]s) 
received at Hanover, N.H., 22 July 1963. 

UT is shown on the curve. The scales are drawn to the origin of the original record. 
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Figure 8. Polar plot of phase and amplitude of NSS (22.3 kc/s) 
received at Hanover, N.H., 20 July 1963. 

UT is shown on the curve. The scales are drawn to the origin of the original record. 



necessary equations and information are given in 
appendix A. 

The use of both of these methods resulted in the 
choice of another origin which is shown in the polar 
plots. The magnitude and direction of the displace- 
ment from the first origin represents the best estimate 
of the phase and amplitude of the groundwave vector, 
or the location of the origin of the skywave vector. 




I 2 3 4 5 6 7 8 9 10" 
UNITS OF CHART DEFLECTION 



Figure 6. Polar plot of phase and amplitude of NSS (22.3 kc/s) 
received at Hanover, N.H., 21 July 1963. 

UT is shown on the curve. The scales are drawn to the origin of the original record. 



during the night, and disturbance effects. Therefore, 
a second method was used in conjunction with the first. 
This method utilizes calculated values of diurnal phase 
and amplitude changes to determine independently 
the area of probability for the skywave origin. The 



4.2. Results 

The linear plots (figs. 9 through 12) were made with 
the new origin and, it is believed, show essentially 
the skywave alone. These contain the same infor- 
mation as the polar plots but in a more convenient 
form. 

The amplitude and phase scales used on the linear 
plots are the same as those used on the original re- 
cords and in the polar plots. Since the amplitude and 
phase values normally vary during the 2 hr eclipse 
period, the midpath noon values averaged over the two 
control days were used as the reference for the cal- 
culations and curves of the eclipse effect. 

From the linear plots, the phase change at any time 
during the eclipse can be determined, and the corre- 
sponding change in path length calculated by (2). 
The height change can then be determined by substitu- 
tion into (1). By this method, it was determined that 
the height change during the eclipse reached a maxi- 
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mum of 10.9 km at the NAA-Hanover one-hop reflec- 
tion point and 6.85 km at the NSS-Hanover one-hop 
reflection point. 

A ratio of amplitudes and a corresponding ratio of 
relative reflection coefficients were then calculated 
from (10) or (11), using the height changes calculated 
above to compute Pn/Pd- The amplitude ratio is 
graphically presented for the NAA-Hanover path 
(fig. 13) and the NSS-Hanover path (fig. 14) for the 
eclipse period. The average values for the two days of 
background data, with the bars on the average values 
giving the range of variation over the two days, are 
shown for comparison. These figures also show, on 
the same time scale, the result of calculations of the 
height change, and the unobscured fraction of the sun's 
visible surface. The daytime skywave and ground- 
wave from NAA are nearly equal in amplitude and 
opposite in phase. The resultant signal is small, 
which causes low amplitude and consequent uncer- 
tainty, since the signal is very sensitive to small iono- 
spheric changes. The NSS signal amplitude and 
phase are such that the time of the NSS amplitude 
maximum is more sensitive to the position of the new 
origin than are the times of the other amplitude and 
phase maxima. 
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Figure 10. Skywave phase and amplitude of NAA (14.7 kc/s) 
recorded at Hanover, N.H., plotted with new origin for 20 July 1963. 
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FIGURE 9. Skywave phase and amplitude of NAA (14.7 kc/s) 
recorded at Hanover, N.H., plotted with new origin for 19 and 
21 July 1963. 



FIGURE 11. Skywave phase and amplitude of NSS (22.3 kc/s) 
recorded at Hanover, N.H., plotted with new origin for 21 and 
22 July 1963. 
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Figure 12. Skywave phase and amplitude of NSS (22.3 kc/s) 
recorded at Hanover, N.H., plotted with new origin for 20 July 1 963. 
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Figure 13. Optical obscuration, calculated height change, and 
reflection coefficient ratio for NAA signal during the eclipse of 
20 July 1963. 

The bars indicate the range of values at these times on the control days. 
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FlGURE 14. Optical obscuration, calculated height change, and 
reflection coefficient ratio for NSS signal during the eclipse of 
20 July 1963. 

The bars indicate the range of values at these times on the control days. 



5. Analysis of Results 

The VLF effects of the eclipse measured at Han- 
over are summarized in table 3. In most cases the 
changes in magnitude and phase of the reflection 
coefficients were not as large as the regular diurnal 
variation. The only possible exception was the change 
in the apparent magnitude of the reflection coeffi- 
cient for the 22.3 kc/s signal. The limited amount of 
data available leaves enough uncertainty about the 
diurnal variation to make the value of the comparison 
of the magnitude of the diurnal and the eclipse varia- 
tions somewhat doubtful. As previously indicated, 
the origin used yields values which agree approxi- 
mately with the VLF propagation calculations. 
However, the magnitude of the deduced diurnal 
variation is sensitive to the position of the new origin. 
Equation (5) indicates very little difference in the 
magnitude of the parallel reflection coefficient, and 
consequently in the amplitude, for night and day. 
The data are consistent with this result. It is difficult 
to draw any conclusions about the meaning of the 
magnitude of these effects because of the uncertainty 
about the behavior of the solar radiation versus time, 
as we shall see later. 

One of the most striking features of the eclipse effect 
is that the maximum signal amplitude occurred con- 
siderably before the optical eclipse maximum. The 
magnitude of these time differences is such that 
they could not be due only to errors in timing. The 
probable explanation for these early maxima will be 
taken up later in the paper. 

The changes in the amplitude and phase of the VLF 
signals during the eclipse are caused by changes in 
the ionization density in the D region. In the case of 
the phase of the signal, we may reasonably assume that 
the signal is reflected from a height at which the ioniza- 
tion density N is approximately 300/cm 3 [Budden, 
1953] . The phase then essentially indicates how the 
height at which this ionization density occurs varies 
as a function of time. The time at which this height 
is a maximum corresponds to the time at which a mini- 



952 



mum in the ionization density is present at a given 
height. 

The time delay of the phase effect relative to the 
maximum optical obscuration was used to calculate 
effective recombination coefficients. These calcula- 
tions are based upon the following assumptions: 
(1) The ionizing radiation which affects the phase was 
a minimum at the time of maximum optical obscuration 
and; (2) the method of analysis given by Appleton 
[1953] may be applied here. Appleton assumes an 
aN 2 type of reaction. This is valid only when the ion 
to electron density ratio X. is constant, i.e., no attach- 
ment type of reaction is present. He shows that for 
the aN 2 type of reaction, when the electron production 
rate has an extreme value, the delay time At for the 
extreme value of N is approximately At=ll2aN y for 
small At, where a is the effective recombination co- 
efficient, and N the electron density. This approxi- 
mation can be better applied to small delay times as 
recorded during the eclipse than to larger ones such 
as during SPA, etc. If k does not remain constant, 
Lerfald et al., [1965] show that the times of minimum 
ionization should not be delayed after the eclipse 
maximum. The fact that a delay is observed may indi- 
cate that X does remain constant or nearly so. In any 
case it means that the change in A does not completely 
dominate the process. 

Using Appleton's expression, and taking the value 
for the electron density at the reflection point of 
300/cm 3 we can calculate the effective recombination 
coefficient for the regions which affect the phase of 
VLF signals given in table 3. The results are 
a = 2.8 X 10" 5 for 14.7 kc/s at a height of about 85 km, 
and a = 1.4 X 10" 5 for 22.3 kc/s at 80 to 82 km (based 
upon the calculated height change and the assumed 
74 km daytime height). 

Now let us consider the time of the maximum ampli- 
tude. The maximum amplitude occurs several min- 
utes before the maximum optical obscuration, 
especially in the case of the NSS measurements; this 
time difference is sufficiently large to rule out timing 
errors as the sole source of the time lead, as mentioned 
previously. 

It appears then that the explanation for the early 
time of the maximum amplitude lies in some physical 
effect caused by the eclipse or in some propagation 
effect. The first explanation that one might consider 
is that of some sort of interference phenomenon. The 



actual ionosphere is, of course, much more compli- 
cated than the single slab model considered so far. 
Studies of ionosphere models with exponential ioniza- 
tion profiles, for instance, have revealed very compli- 
cated behavior of the reflection coefficient as a function 
of changes in the profile [Wait and Walters, 1963] . It 
is possible that such a configuration might produce, 
merely by interference between partially reflected rays, 
an interference pattern which has a maximum ampli- 
tude at any time. The fact that the maximum ampli- 
tude occurred considerably before the optical 
maximum of the eclipse on both VLF paths makes it 
seem unlikely that this could be merely a fortuitous 
interference effect. When the times of occurrence 
of the maximum amplitude are analyzed further, we 
shall see that the particular times at which the maxi- 
mum occurred have a significance which is unlikely 
to have occurred by mere accident. 

Let us then consider another possible explanation 
for this early time of maximum amplitude. Assume 
that the ionosphere may be represented, as far as the 
effect is concerned, by two homogeneous slabs. The 
lower slab, of arbitrary thickness, has values of elec- 
tron density and collisional frequency such that it 
produces only attenuation of a wave passing through 
it. The upper, semi-infinite slab has the proper elec- 
tron density and collisional frequency so that it 
produces a reflection at the lower edge. It is possible 
that the early time of maximum amplitude then could 
be caused by some component of the ionizing radiation 
which affected only those lower heights where the 
majority of the attenuation, and thus the change in 
the amplitude or the apparent magnitude of the re- 
flection coefficient, is produced. 

The variation in the amplitude is caused by the 
variation in the ionization density present in part of 
the D region where the attenuation is produced. The 
time variation of the amplitude then is indicative of the 
time variation of the ionization density and the delay 
time may be treated in a similar manner. 

The question of the source of the ionizing radiation 
immediately comes to mind. Bracewell [1952] has 
stated that he saw no indication of any effect of the 
obscuration of active areas on the VLF signals during 
the eclipse of 28 April 1949. No other mention of 
this effect in the D region has been found in the litera- 
ture, although it has been reported in ^-region studies- 
[Minnis, 1955; Ratcliffe, 1956]. 



Table 3. Summary of VLF effects of 20 July 1963 eclipse measured at Hanover, N.H. 







Solar obscuration 


Maximum 


Station, location, frequency, and distance 


Magnitude of reflection 
coefficient ratio 


Change of phase and apparent reflection 
height 




UT 


% Mag 


UT 


Ratio 


% of Diurnal 
variation 


UT 


Phase 


Height 


% of Diurnal 
variation 


NAA, Cutler, Me. 14.7 kc/s, 406 km 
NSS, Annapolis, Md. 22.3 kc/s, 633 km 


2145 
2149 


100 
88 


2143 
2142 


1.54 
1.83 


77 
100 


2146 
2151 


Degrees 
144 

102 


km 
10.9 

6.85 


68 
43 
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The ionizing radiation which produces the daytime 
D region is generally considered to consist of Lyman 
Alpha (1217 A), x rays of X < 10 A [Nicolet and Aiken, 
1960] and cosmic rays [Reid, 1964]. The cosmic 
radiation should provide a constant ionizing source 
which will not vary during the eclipse. It appears 
that at least as a first assumption, cosmic rays will not 
contribute any asymmetry about the time of maximum 
eclipse so they will not be considered further. Fried- 
man [1963] found from rocket measurements during 
the 12 October 1958 eclipse that the Lyman alpha 
radiation fell to values comparable to the general sky 
background level, but about 12 percent of the x-ray 
flux persisted throughout totality. Furthermore, 
Friedman found that the x-ray intensity in the corona 
was six times as great on the west limb, which con- 
tained several active plages, as on the east limb, which 
was free of activity. 

Another measurement by Friedman [1963] with a 
rocket-borne pinhole camera in April 1960 showed that 
the x-ray emission extended about 0.06 solar radii, or 
43,000 km, above the solar limb. The examination of 
available x-ray records and Calcium-K line photos led 
Friedman to the conclusion that most x radiation comes 
from the corona above bright plages. A further com- 
parison shows that the x-ray source contours very 
closely resemble the contours of the 9.1 cm spectro- 
heliograms made at Stanford University [CRPL 
F-Series 1963]. 

The optical eclipse is caused by the obscuring of the 
sun's visible disk by the disk of the moon. These 
two disks are approximately the same size. It is 
evident then that an area above the visible solar disk, 
located at or near the limb of the sun, would extend 
considerably beyond the visible disk. If it is near the 
solar equator, it can be eclipsed and uncovered earlier 
or later than the first or last contact of the optical 
eclipse, depending upon whether this region was 
located near the east or west limb of the sun. These 
considerations suggest the possibility that the eclipsing 
of active x-ray producing areas above the optical disk 
should be considered as an explanation for the early 
time of maximum amplitude. 

A map of the sun's surface at the time of the eclipse 
[Billings et al., 1963] shows that the visible disk covers 
312 to 134° heliographic longitude. There was an 
active 5303 A emission region located between about 
130 and 165°, at latitudes between about 3° S and 12° N. 
Two calcium plages were noted at 139.5° longitude and 
14° N latitude, and 149.5° longitude and 4° N latitude. 
This region also appears on the 9.1 cm spectrohelio- 
gram made at Stanford University at the time of the 
eclipse. The region would have been at the west limb 
of the sun. The active x-ray producing area above 
this region would have been eclipsed before first con- 
tact of the optical eclipse, and would have been 
uncovered before optical maximum. 

A rough calculation was made, assuming that (1) the 
apparent size of the solar and lunar disks are the same, 
(2) the rate of motion was uniform, and required 2 hr 
from first to last contact, and (3) the upper limit of 



0.06 solar radii found by Friedman is typical for the 
maximum height of x-ray sources above the limb. 
The calculations indicated that the region above the 
plages and active regions should be uncovered 3.7 
min and 7.2 min before optical maximum for the NAA 
and NSS one-hop reflection points, respectively. Let 
us consider the NAA signal first. The maximum 
amplitude occurred 2 min before the optical maximum. 
Therefore, the maximum amplitude occurred about 
1.7 min after the calculated time of the uncovering 
of this active area. The maximum amplitude of the 
NSS signal occurred about 7 min before optical maxi- 
mum. Therefore, the maximum amplitude occurred 
about 0.2 min after the calculated time of the uncover- 
ing of the active area. The corresponding delay times 
for the maximum phase effect on NAA and NSS, 
respectively, were 1 and 2 min. These delay times 
all he close together, with the calculated delay for the 
amplitude effect on the NAA signal lying in between 
the two values calculated for the phase effect. This 
suggests that at least the amplitude delay for the NAA 
signal is fairly reasonable. As was mentioned before, 
the time of the NSS signal amplitude maximum is more 
dependent than the others on the position of the new 
origin for the polar plot. The fact that the measured 
time delay is only 0.2 min is believed to be largely due 
to this fact. In addition, the signal-to-noise ratio of 
the NSS signal was a bit smaller than that of the NAA 
signal which contributes a little more uncertainty to 
the time of maximum. 

It is very difficult to determine the magnitudes of 
the errors' involved in the measurement of these time 
delays. Actual timing errors on the record are ex- 
pected to amount to not more than ±0.5 min. The 
amplitude and phase measured on the record, however, 
is the result of an effect which is integrated over the 
area of the first Fresnel zone. The time of optical 
maximum varies over this first Fresnel zone by about 2 
min in the case of the NAA signal and by about double 
this in the case of the NSS signal. Since the actual 
contribution of any part of these areas decreases as 
its distance from the center increases, it is very difficult 
to set an actual numerical limit on the effects of this 
variation on the time of maxima. Although the delay 
times measured on the 22.3 kc/s signal from NSS tend 
to be a bit higher than those measured for the 14.7 
kc/s signal from NAA, there is not sufficient precision 
available to say that the difference is significant. 
It is better to find a reasonable average value for the 
results. As indicated in table 4 this yields a value of 
approximately 2 X 10~ 5 for the effective recombination 
coefficient. It seems very unlikely that the error of 
the measurement could account for the difference of 
roughly a factor of ten from many previous values 
obtained by other observers. 

The values determined by others for the effective 
recombination coefficient are also tabulated in table 
4 [Mitra and Jones, 1954]. The other values were not 
taken during an eclipse but during daytime conditions 
or during SPA's when X would be expected to be 
constant. It should be noted that the delay times 
found during most of these are quite long, ranging all 
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the way from 4 § to 80 min. Bracewell [1952] deter- 
mined that during an earlier eclipse the 'delay time 
could be not larger than 6 min. In addition, Bracewell 
and Straker [1949] stated that the relaxation time 
determined from SPA's due to solar flares does not 
exceed 5 min and may be much less. They further 
state that in at least one case the time is 2.5 min or less. 
Furthermore, Gallet [1955] found that absorption 
measurements made during the eclipse of June 20, 
1955, yield a time delay of not more than 1 or 2 min. 
These time delays are very close to the eclipse values 
determined in this paper and would yield similar effec- 
tive recombination coefficients. 

There is therefore considerable experimental evi- 
dence of these high effective recombination coeffi- 
cients. The values given by Mitra and Jones [1954] 
bracket the daytime values (=1X10~ 7 at 80 km) 
shown by Crain. The values calculated here are 
nearer to the values (= 2 X 10" 6 ) given by Crain for the 
nighttime case. It is also possible that the ionospheric 
conditions in the D region at the eclipse maximum 
were close to those of nighttime, since the apparent 
reaction time is about the same as the duration of 



totality (2 min). These results would be reasonable 
for the ionosphere about 20 km lower, i.e., 60 km 
[Mitra and Jones, 1954; Crain, 1961]. It is perhaps 
more likely that this indicates an attachment type of 
reaction is actually present, although not dominant. 

The reversal of the direction of change of the height 
and reflection coefficient ratio for the NAA signal and 
the change in slope of the reflection coefficient ratio 
for the NSS signal, which occurred several minutes 
after the optical maximum, could have resulted from 
the obscuration of some active region. However, 
these times do not appear to correspond to the covering 
or uncovering of any obvious feature on the sun. 
It appears more likely that it is due to a relatively weak 
interfering signal such as the second order ionospheric 
reflection. 

It seems a little strange at first that there was no 
noticeable effect on the signals when the x-ray produc- 
ing area was first obscured a few minutes before optical 
first contact. The answer probably lies in the fact 
that the fractional change in radiation was small, 
since the visible disk and any other x-ray producing 
areas were still unobscured. 



Table 4. Comparison of calculated values of the effective recombination coefficient 
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1.4- 10" 
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Mitra and Jones 


Benner 


80 
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150 


60 


1.4- 10" 
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80 


Polarization 
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80 
80 


1 • 10" 4 
1 • 10" 
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Mitra and Jones 


Nearhoof 


90 


Phase height 


150 


30 


5.2 • 10" 


3060 
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Mitra and Jones 


Jones 


75 


"D«" phase ht. 


16 


25±5 


3.6 • 10" 


280 
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Bracewell 
and Bain 
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1.9- 10- 6 


Mitra and Jones 


Ellison 


70 
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2.2 • 10- 6 


Mitra and Jones 


Ellison 


73 


SPA decay 


16 








6.3 • 10-« 


Mitra and Jones 


Bracewell 


65-75 


SPA decay Cls 
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16 


8 






3.5 • 10" 6 


Mitra and Jones 


Mitra and Jones 


65-75 


SPA decay Cls 
3 — 3 + flares 
mean values 


16 


4.5 






6.3 • lO" 6 


Mitra and Jones 


Mitra and Jones 


80-85 


Amplitude and 


14.7 




[ ~ 6 • 10~ 3 


-300 


- 2 • 10- r> 


Crary and 


Crary and 




phase during 


22.3 


2 


Schneible 


Schneible 




eclipse 

















6. Conclusions 

The measuring techniques and the methods used to 
analyze the VLF records from the eclipse yield re- 
sults which are consistent with the theories of propaga- 
tion and the present knowledge of ionospheric 
conditions during the eclipse. The effective recom- 
bination coefficient calculated from these eclipse 
data is higher by a factor of 10 or more than the values 
calculated by most experimenters. In addition it is 
nearer the value calculated from theoretical considera- 
tions by Crain for the nighttime D region than it is 
to that calculated for the daytime D region. Other 
observers have obtained values of time delay from 
eclipse and SPA data which are of the same magnitude 
as those obtained in this paper. It is possible that 



the actual heights involved are considerably lower, 
where the recombination coefficient is greater, but it 
is more likely that this indicates an attachment type of 
reaction rather than true recombination [Reid, 1964]. 

The fact that the NAA amplitude delay time agrees 
so nicely with the delay times calculated from the 
phase results indicates that the assumption about the 
x-ray source is very likely the true explanation. In 
addition these assumptions are in line with the known 
facts about the Lyman alpha, cosmic rays and x radia- 
tion which are the source of Z)-region ionization. 

It is plausible then, that the ionization at the higher 
heights which affects the phase (say, above 80 km) did 
result from the Lyman alpha radiation from the visible 
disk, while the ionization at the lower heights, which 
affects the amplitude, was caused by the x radiation. 
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This radiation can reasonably be assumed to arise from 
the area above the visible disk over an active area 
at the west limb. This active area was also observed 
in 5303 A emission [Billings et al., 1963] and also in 
the 9.1 centimeter spectroheliograms [CRPL-F 
Series, 1963] made of the sun around the time of the 
eclipse. 
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7. Appendix A. Calculation of the Ex- 
pected Diurnal Changes in Phase and 
Amplitude 

The distance of a single trip between earth and iono- 
sphere is given [Helliwell, 1953] by 

P 2 = 2R(h' + R)l[l-cos(180DlTrR) + h' 2 l (1) 

where P = ray path distance in km, 
h' = height in km, 
D = total ground distance in km, 
R = mean earth radius in km. 

Then, the expected diurnal phase change is found by 
substituting the values for the night and day ray path 
distances into the expression 



A0 = 



2(P N -P D )(360°) 



(2) 



where the N and D refer to night and day conditions, 
respectively. The values of 74 km for the noon reflec- 
tion height and 16 km for the diurnal change have 
been indicated by others [Bracewell et al., 1951] 
and were used here. The change in height for a given 
phase change is relatively insensitive to the value of 
height assumed in this range, so the exact value of the 
height is not critical. The assumed reference height 
of 74 km is used only to define the height range of the 
effects being measured. 

Determining the theoretical diurnal amplitude 
change requires two simplifying assumptions. First, 
it is assumed that the ionosphere is a sharply bounded 



homogeneous medium with a constant magnetic field; 
second, that the Q—L approximation is valid [Crary, 
1961]. Budden [1951], and Wait and Murphy [1956] 
give expressions for the refractive indices /Xo, \ie, and 
reflection coefficients of such a medium in terms of the 
parameters co r and r. 

Reasonable values of these parameters for the 
ionosphere for VLF have been determined from experi- 
mental data by Wait and Murphy [1956, 1957], Budden 
[1951], Bracewell et al. [1951], and Crary [1961]. 
These values are: for night a) r = 5 X 10 5 , r = 60°, and 
A' = 90 km; for day w r = 1.5Xl0 5 , t=15°, and h! 
= 74 km. Then the angle of incidence at the iono- 
sphere, as expressed by Helliwell [1953], is 



= tair 



{-C3VK- 



/180DY 

cos l^FJ 



(3) 



Budden [1951] gives equations for the reflection 
coefficients in the form 



\\Rl =2i C()LtoC — fl e C e )/D 



(4) 



\\R\\= [(Mo + lie) (C 2 — C Ce) 

+ ( M oiLte-l)(Co+Ce)C]/Z> (5) 

where \x e sin d e = fiio sin 6 = fi sin 6 

D = (fJL + fie) (C 2 + CoCe) + (We + 1) (Co + C e )C 

C = cos 6 

Co = cos O 

C e = COS e . 

If we assume that the transmitting antenna has a 
cosinusoidal pattern in elevation over a perfectly 
reflecting ground, the radiation field is given by 

E\\=KVP cos ft, (6) 

where E\\ = the normal or parallel field component, 
K=sl constant of proportionality, 
P = the radiated power in kw, 
ft = the elevation angle at the ground. 

The response of a loop antenna over a highly con- 
ducting ground in terms of the electric field [Crary, 
1961] is 



v i ~ h' L (E\ | cos cj> — E± sin </> sin ft), 



(7) 



where i;l = the induced voltage in the loop, 

h\ = the effective height of the loop over a highly 

conducting ground, 
<f> = the angle between the plane of the loop and 
the plane of incidence. 



The abnormal loop response is then given by 
KVF cos p ift L sin ft, 



v A =- 



Pn 



(8) 
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and the normal loop response is given by 



v N = - 



K VF cos ft ||fi|| 

Pn 



(9) 



The night to day amplitude ratios, A\\nIA\\D and 
A±nIA±d, are then found from (8) and (9) to be 



and 



(10) 



(11) 



A\\ N = cos Pn\\R\\n Pd 
A\\ D cos Pd\\R\\d Pn 



A± N _ sin Pn cos @n\\R]_n Pd 
A^ D sin fi D cos (3 d \\R]_d Pn 
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